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ABSTRACT 

Context. We investigate the stellar population and the origin of diffuse light around brightest cluster galaxies. 

Aims. We study the stellar population of the dynamically hot stellar halo of NGC 331 1, the brightest galaxy in the Hydra I cluster, 
and that of photometric substructures in the diffuse light to constrain the origin of these components. 

Methods. We analyze absorption lines in medium-resolution, long-slit spectra in the wavelength range 4800 - 5800 A obtained with 
FORS2 at the Very Large Telescope. We measure the equivalent width of Lick indices out to 20 kpc from the center of NGC 331 1 
and fit them with stellar population models that account for the [a/Fe] overabundance. 

Results. Stars in the dynamically hot halo of NGC 3311 are old (age > 13 Gyr), metal-poor ([Z/H] 0.35), and alpha-enhanced 

([ff/Fe] ~ 0.48). Together with the high velocity dispersion, these measurements indicate that the stars in the halo were accreted from 
the outskirts of other early-type galaxies, with a possible contribution from dwarf galaxies. We identify a region in the halo of NGC 

331 1 associated with a photometric substructure where the stellar population is even more metal-poor ([Z/H] 0.73). In this region, 

our measurements are consistent with a composite stellar population superposed along the line of sight, consisting of stars from the 
dynamically hot halo of NGC 3311 and stars stripped from dwarf galaxies. The latter component contributes < 28% to the local 
surface brightness. 

Conclusions. The build-up of diffuse light around NGC 3311 is on-going. Based on the observed stellar population properties, the 
dominant part of these stars may have come from the outskirts of bright early-type galaxies, while stars from stripped dwarf galaxies 
are presently being added. 



Key words, galaxies: cluster: individual: Hydra I (Abell 1060) - 
galaxies: halos - galaxies: formation 
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1. Introduction 

Brightest cluster galaxies (BCGs) are the most massive and 
luminous objects located at the center of galaxy clusters. 
According to the hierarchical formation scenario, they assemble 
through the merger of smaller objects that are dragged into the 
cluste r's center by dynamical friction (e.g. iDe Lucia & Blaizol 
l2007h . During the evolution of the cluster and the formation 
of the central BCG, stars are stripped from their host galax- 
ies. Part of the stripped material is captured by the halos of 
the BCG and other massive galaxies, while the other part re- 
mains free-floating in the cluster potential, forming the so-called 
intra-cluster light (ICL). The details of the interplay between the 
formation of the BCG halos and the ICL are not yet fully un- 
derstood. Observations of both components can give key infor- 
mation about the processes involved in the evolution of cluster 
galaxies. 

In the last few years, progress was achieved in the char- 
acterization of the stellar population properties of the halos 
of BCGs at galactocentric radii larger t han three effective 
radii by means of lo ng-sIit spectroscopy jBrough et al.l 120071 
ICoc cato et al. 2010a b) or photometric studies of globular clus- 
ters (Forbes et al. 2011; Arnold et al. 2011). Stai's in the BCG 
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halos are generally old (>10 Gyr), they have solar or sub- 
solar metallicity and super-solar a-enhancement. Sometimes a 
change in the radial profiles of metallicity, c-enhancement or 
color is observed, indicating that the stellar population proper- 
ties in the outskirts are different from those in the central re- 
gions. The emerging interpretation is that the inner and outer re- 
gions formed in two separate phases: the central parts formed 
"in situ" (or through a single major merger), while the halo 
was accreted later through a series of minor mergers. This 
"two stages" formation process for BCGs is also supported 
by numerical simulations (e.g. Abadi et al. 2006; Oser et aT| 
2010) and photometric studies that compare the BCGs scal- 
ing relations between color gradients, velocity dispersions, and 
sizes with thos e of other early-type galaxies of similar mass 
or luminosity (iBernardi et al.l I2OIII) . Minor merger accretion 
in the second stage of this scenario is also consistent with 
the observed size evolution of early-type gala xies over redshift 
(iTruiillo et al.ll2007t Ivan Dokkum et al."2010b. hydrodynamical 
cosmol ogical simulatio ns (Naab et al. 2009), and spectroscopic 
studies (ICoccato et al.ll2010bt ITruiillo et aI.II201 11) . 

It is not clear yet whether the major contribution for the 
formation of the ICL comes from stars stripped in galaxy 
groups and subsequently captured by the clus ter potential 
(Rudicketal. 2006, 2009; Kapfe reret al.1 l2010h . from stars 
of galaxies th at interact w i th each other and with th e clus- 
ter tidal field dMoore et al.1 Il996t IWillman et all |2004 . from 
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Fig. 1. Digital Sky Survey image of the central regions of the 
Hydra I cluster The two brightest galaxies NGC 3311 and 
NGC 3309 are labeled. The thick black line represents the loca- 
tion of the FORS2 long slit. The orange circle marks the slit cen- 
ter on top of the dwarf galaxy HCC 26. The spectrum of HCC 26 
was extracted from the central 6" (see text for details). The col- 
ored lines identify the location of the five slit sections (bins A-E) 
where we extracted the spectra of the halo of NGC3311. Thin 
continuous, dashed and dotted lines represent a pproximately the 
contou rs of a photometric substructure found bv lArnaboldi et al.] 
(1201 Ih . We adopt a distance to NGC 331 1 of 51 Mpc, from the 
NASA/IPAC extragalactic database. 



stars stripped from galaxies merging into the central BCG 
(Murante et al. 2007; Puchwein et al. 2010), or fr om entirely dis - 
rupted low surface brightness galaxies or dwarfs (lGnedinll2003l) . 

The connection between the accretion of the BCG halo and 
the formation of the ICL is hard to probe from an observational 
point of view, given th e low surface brightness of these regions. 
IWiUiams e t ail(l2007h measured the stellar population parame- 
ters for a sample of ~5300 intra-cluster stars in one field in the 
Virgo cluster Although they found a heterogeneous stellar popu- 
lation, 70%-80% of the observed intra-cluster stars are old (age 

> 10 Gyr) and metal-poor (mean [M/H] 1.0). On the other 

hand, other studies report super-solar me tallicities in the di ffuse 
light in Abell 3888 (I.OZq < Z < 2.5Z,.. lKrick et al.ll2006h that 
are interpreted as ICL formed early with the collapse of the main 
cluster, unrelated to the BCG halos. From a kinematic point of 
view, the line-of-sight velocity distribution of intra-cluster plan- 
etary nebulae in one Coma cluster core field is consistent with 
the ICL being formed from stars stripped from the halos of the 
two brightest galaxy in the cluster while they orbit each other 
dGerhard et al. 2007), while recent studies of absorption line 
spectra in the RX J0054. 0-2823 cluster are consistent with the 
ICL being the remnant of tidally destroyed galaxies and stream- 
ing in the central regions of the cluster, in w hich the three centra l 
giant ellipticals act as "grinding machines" (iToledo et al.ll201 11) . 
The origin of the ICL and its connection with the BCG halos 
is thus not fully understood, and it is likely that several mecha- 
nisms are acting in its formation. 



In this context we have been carrying out a project aimed at 
studying the formation of the intra-cluster stellar component in 
the Hydra I galaxy cluster (Abell 1060), and this paper is the 
fourth in this series. We observed a rapid increase of the ve- 
locity dispersion radial profile in the ha lo of NGC 3311, the 
BCG in Hydra I dVentimi^ha et al.ll201()l) . This was interpreted 
as caused by a dynamically hot stellar halo around NGC 3311, 
dominated by intra-cluster stars. We also observed the veloci- 
ties of planetary nebulae in the core of the Hydra I cluster, the 
major component of which is consist ent with bei ng part of this 
halo of intra-cluster stars (Ventimiglia et al.ll201 lb . The distribu- 
tion and kinematics of the entire sample of intra-cluster plan- 
etary nebulae indicat e that most of the diff'u se light is still not 
phase-mixed. Finally, lArnaboldi et al.l (1201 Ih found that the hot 
stellar halo is not centered on NGC 3311, and detected a pho- 
tometric sub-structure (excess of light, hereafter) superposed on 
the intracluster halo. The contours of light that best approximate 
this photometric excess are shown in Figure [T| Their results in- 
dicate that this excess of light has line-of-sight velocity ~ 1200 
km s ' relative to the stellar halo around NGC 331 1, and that it 
is related to one or more dwarf galaxies disrupted in this region. 
The build-up of the ICL and the dynamically hot stellar halo 
around NGC 3311 are therefore ongoing through the accretion 
of material from galaxies falling into the cluster core and tidally 
interacting with its potential well. 

In this paper we investigate the nature of the dynamically 
hot stellar halo around NGC 3311 and the excess of light by 
studying their stellar populations. Section |2] presents the spec- 
troscopic observations and the data reduction steps. Section |3] 
describes measurements of equivalent widths of the Lick indices 
and the derivation of stellar age, metallicity, and or-enhancement. 
In Section |4] we discuss the origins of the stellar halo of NGC 
3311 and the excess of light, and finally in Section|5]we conclude 
with a brief summary. 

2. Spectroscopic observations and data reduction 

The long-slit spectra were obtained during the nights of March 
25-28-2009, wi th the focal reducer and l ow-dispersion spectro- 
graph (FORS2, lAppenzeller et al. I [19981) mounted on the Unit 
Telescope 1 of the European Southern Observatory (ESO) in 
Paranal (Chile). The long-slit was oriented at position angle 
PA - 142°, and centered on the dwarf galaxy HCC 26, which 
is ~ 45" northeast from NGC 3311 (Fig.[B- The instrumental 
setup had the I'.' 6 wide long-slit and the 1400V+18 grism with 
instrumental spectral resolution of FWHM ~ 210 km s"'. The 
total exposure time wa s four hours. For m ore information on the 
instrumental set-up see I Ventimiglia et al ](|2010). 

Data reduction (overscan and bias subtraction, flat-fielding 
correction, wavelength and flux calibrations) was performed 
with the FORS2 reduction pipeline version 3.0 under the ESOrex 
environment and IRAlQ and IDL ad hoc scripts. The contri- 
bution from scattered light in the instrument was interpolated 
from the portions of the F0RS2 detectors not illuminated by 
the long-slit (typically a few rows on the CCDs borders) and 
then subtracted. Scattered light amounted to 7 + 3 % of the sky- 
background value, with negligible gradients along the dispersion 
direction. Wavelength calibration and spectral rectification were 
performed using He, Hg+Cd, and Ne calibration lamp spectra 

' Updated version of the F0RS2 pipeline and ESOrex can be found 
at http://www.eso.org/sci/software/pipelines/ 

- IRAF is distributed by NOAO, which is operated by AURA Inc., 
under contract with the National Science Foundation. 
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Fig. 2. Spectra extracted in the selected slit sections (bins A-E) and in the slit center (HCC 26). Black: observed spectrum, red: 
best-fit model, green wavelength region excluded from the fit. 



and were cross-checked with the most intense emission lines in 
the sky spectrum. Errors on the wavelength calibration were ~ 5 
km s Cosmic rays were removed using the FILTER/COSMIC 
routine in MIDA^. Long-slit spectra on blank fields in the Hydra 
I cluster were also obtained close in time to the scientific spec- 
tra to evaluate the sky background and to correct for large-scale 
illumination patterns from the non-uniform illumination along 
the slit. All spectra were corrected for atmospheric extinction 
and flux-calibrated using spectro-photometric standard stars ob- 
served during twilight as part of the FORS2 standard calibration 
plan. 

2.1. Sky subtraction 

Accurate sky subtraction is important when measuring spectra 
of low surface brightness regions. To limit systematic effects as 
much as possible, we accounted for the sky-subtraction in two 
ways, following the strategy adopted in C occato et al.. C2010a) . 

' MIDAS is developed and maintained by ESO. 



The mean flux of the continuum of the sky spectrum extracted 
from the outermost portions of the long slit was ~ 25% higher 
than that of the sky spectrum extracted from blank fields in the 
Hydra I cluster. This indicated that the contributions from the 
halo of NGC 3311 and from intra-cluster stars were not negli- 
gible in the long-slit ends. We therefore decided to use the sky 
spectrum obtained from blank fields. Although they were taken 
close in time to the scientific spectra, they were not simultane- 
ous. To minimize the effects of the variation in time of the in- 
tensity of the sky em ission lines, we proceeded in a way sim- 
ilar to that of IDaviesI (12007). For each night's sky spectrum, a 
fourth order polynomial was fitted to the continuum and then 
subtracted. The mean flux of the sky-background continuum in 
the 4800 - 5600 A range did not vary by more than +5% from 
night to night. The residual spectrum that contained emission 
and absorption line features was multiplied by a constant factor 
and added to the sky continuum previously subtracted. The con- 
stant factor (typically on the order of 1.20 - 1.30) was chosen to 
minimize the difference between the intensities of some selected 
sky emission lines measured in the scientific spectra and those 
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Table 1. Line strengths of Lick indices in the hot stellar halo around NGC 33 11. 



ID 


D 


PA 


H/3 


Fe5015 




Fe5270 


Fe5335 


Fe5406 


Fe5709 




["] 


n 


[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


(i) 




(3) 


(4) 


\->) 


(OJ 


U) 


(o) 


(9) 




D 


42.8 


47.1 


0.77 ± 0.09 


3.14 ±0.19 


3.02 ±0.12 


1.63 ±0.10 


1.84 ±0.17 


1.13 ±0.09 


0.81 ±0.08 


C 


50.2 


83.7 


0.96 ± 0.09 


3.42 ±0.18 


3.41 ±0.10 


1.69 ±0.10 


2.05 ±0.15 


1.26 ±0.09 


0.85 ± 0.08 


E 


57.1 


10.3 


0.66 ±0.10 


2.57 ± 0.22 


2.57 ±0.18 


1.24 ±0.11 


1.58 ±0.21 


0.77 ±0.12 


0.75 ±0.10 


B 


67.2 


102.5 


1.11 ±0.09 


3.40 ±0.19 


3.18±0.11 


1.78 ±0.10 


1.80 ±0.15 


1.23 ± 0.09 


0.81 ±0.08 


A 


108.4 


118.8 


0.99 ±0.10 


4.05 ± 0.24 


3.06 ± 0.26 


1.68 ±0.18 


1.61 ±0.23 


0.77 ±0.13 


0.69 ±0.11 


HCC26 


45.0 


64.0 


1.85 ±0.09 


2.53 ±0.18 


1.85 ± 0.07 


2.05 ± 0.07 


1.87 ±0.13 


1.22 ± 0.07 


0.61 ± 0.09 



Notes - Col. 1: Slit bin name, according to Fig.[T] Col. 2: Distance of the bin center from the center of NGC 331 1 Col. 3: Position angle of the 
bin center with respect to the center of NGC 331 1. North is 0°, East is 90°. Cols.4-1 1: Equivalent widths of the Lick indices and their errors. 



measured in the sky spectra. The entire process was optimized 
in the wavelength range 4800 A - 5800 A. 

Finally, the sky-subtracted spectra were coadded using the 
center of HCC 26 as reference for alignment. 

HCC 26. The spectrum of HCC 26 was extracted from the 
central 6" of the FORS2 long-slit (orange circle in Fig. [T]i. The 
contribution of NGC 3311 was interpolated from the adjacent 
~5'.'8, where residual flux of HCC26 is negligible (lower than 
4% of the HCC 26 central value and comparable to the local 
noise level) and removed from the spectrum. 

NGC 3311 halo. The spectrum of the halo of NGC 33 1 1 was 
extracted along five sections of the FORS2 long-slit (bins A-E) 
from the sky-subtracted spectrum. Their locations are displayed 
in Fig.[T] The edges of the central bins C and D are ~ 6" away 
from the slit center, at positions where the contamination from 
HCC 26 is lower than the background noise level. The sizes of 
the bins are chosen to ensure a minimal signal-to-noise ratio ~ 
25 per spectral pixel in each bin, and to account for the CCD 
bad cosmetics and faint foreground stars in the slit. Outside the 
explored region the surface brightness level is too small to permit 
reliable measurements. 

We quantified the residual contamination from sky emission 
lines on the final co-added science spectra to be lower than 2% 
of the subtracted sky itself, and lower than the noise level. This 
was determined in the following way: 

1 . We subtracted the best-fit models from the observed galaxy 
spectra, obtaining the residual spectra R{X) for each slit bin. 
We computed the average noise < rms > values as standard 
deviations of the R{A) spectra in the wavelength range 4900 
-5550 A. 

2. We cross-correlated the R{A) with the continuum-subtracted 
sky spectrum 5 ems (/I) to detect evidence of contamination 
from residual emission lines. No significant correlation was 
found, meaning that the contamination was below the noise 
level. 

3. We constructed a noise spectrum N{X) with mean and stan- 
dard deviation equal to the maximum < rms >, which is 
measured at the outermost slit section (bin A, Fig.|2]). 

4. We determined the minimum value of K that ensured a sig- 
nificant correlation between the 5ems('t) spectra and N{A) + 
K ■ 5em,s('t)- We found K - 0.02, meaning that the high- 
est possible contamination from residual sky-emission lines 
is 2% of the sky value. Indeed, if the contamination were 
higher, it would have been detected in step 2. The 2% upper 
limit will be used in Section [TTI to quantify the maximum 
systematic eiTors caused by residual sky contamination in 
the derived stellar population parameters. 



3. Line strength indices and single stellar 
population 

The spectra extracted from the slit sections were fit- 
ted with the MILES librar y of stellar template spectra 
jSan chez-Blazquez et al.' 2006) using the penalized pixel-fitting 
routine by Cappellari & Emsellem (2004) to recover the stellar 
velocity (V), velocity dispersion (cr), third and fourth Gauss- 
Hermite moments (/i3 and h4) of the line-of-sight velocity dis- 
tribution (LOSVD). In the fitting process, the spectra of the 
MILES library were convolved with a Gauss-Hermite function 
to match the instrumental FORS2 spectral resolution. The ap- 
propriate function was determined by fitting Gauss-Hermite mo- 
ments to the correlation function between the spectra of one 
star observed during the run and its corresponding spectrum in 
the MILES library. The use of a Gauss-Hermite function in- 
stead of the commonly used Gaussian function ensured more 
reliable measurements, because the FORS2 instrumental line- 
spread function slightly deviates from a Gaussian profile. The 
extracted spectra are shown in Figur e[2ltogether with their b est- 
fit models. We refer the reader to IVentimiglia et al.l (1201 Oh for 
the discussion of the kinematics. In this paper we focus on the 
measurements of the Lick indices and the derivation of the stellar 
population properties. 

Galaxy spectra were set to the rest wavelength using 
the measured radial velocities and then broadened to match 
the Lick instrumental spectral resolution (8.4 A at 5100 A, 
IWorthev"& Ottaviani 1997) . The equivalent widths of the line 
strengths of s everal Lick ind i ces we re measured according to the 
definitions of IWorthev et al.l (Il994l) . and their er rors were evalu- 
ated using the empirical formulas bv lCardiel et a l. ( 1998). 

The measurements were also corrected for kinematic broad- 
ening caused by the line-of-sight velocity distribution, using 
the multiplicative coefficients tabulated in Kuntschner (2004!), 
which account for cr, /i3 and /i4. Because no coefficient is tab- 
ulated for the Fe5709, the corresponding factor was determined 
as the ratio of i) the intrinsic values measured on the optimal 
stellar template; and ii) the values measured on the optimal stel- 
lar template convolved with the measured line of sight velocity 
distribution. 

A set of seven Lick standard stars were also observed with 
the same instrumental configuration. The line strengths of the 
Lick indices were meas ured for these s tars an d compared with 
the values tabulated in IWorthev et all d 19941) to measure off- 
sets to the Lick system. Small offsets were observed for the H/3 
(-0.11 + 0.08 A) and Mg b (-0.13 + 0.07 A) indices and were 
applied to our measurements. 
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Table 2. Luminosity-weighted single stellar population parame- 
ters in the hot stellar halo around NGC 33 11. 



ID 

(1) 


a 

["] 

(2) 


Age 
[Gyr] 
(3) 


[Z/H] 
(4) 


[a/Fe] 
(5) 


D 
C 
E 
B 
A 


43.2 
50.8 
61.3 
70.2 
117.0 


13 5+"-" 

13 5+''-"' 

13 5+°-" 
1 J.J „ „ 

13 5+"-" 

13 5+"-" 


-0.46 ± 0.04 
-0.33 + 0.03 
-0.73 + 0.06 
-0.39 ± 0.04 
-0.34 ± 0.05 


0.45 ± 0.04 
0.48 ± 0.03 
0.50 + 0.06 
0.44 ± 0.04 
0.50 ± 0.03 


HCC 26 


45.4 


13 5+"-" 

1 J.J j Q 


-0.85 ± 0.03 


-0.03 ± 0.05 



Notes - Col. 1: bin name, according to Fig.[Il Col. 2: distance of the bin 
center from the center of NGC 3311 projected onto the galaxy major 
axis assuming an axial ratio fo/a = 0.89 (Arnaboldi et al. 2011). Cols. 
3-4: single stellar population best-fit values and their errors. [Z/H] and 
[cf/Fe] are in logarithm of solar units. 



Table [T] lists the final equivalent widths of the Lick indices 
we measured, corrected for kinematic broadening and offset to 
the Lick system. Errors in the equivalent widths also include the 
uncertainties in the offsets to the Lick system and in the correc- 
tion for kinematic broadening. 
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Fig. 3. Equivalent luminosity-weighted single stellar population 
values of age (upper panel), metallicity (central panel) and a- 
enhancement (lower panel) obtained in the five long-slit sec- 
tions, as function of the distance to NGC 3311. Labels in the 
upper panel identify each radial bin. The shaded yellow regions 
indicate the range of SSP parameters that would result if system- 
atic residual effects in the sky subtraction at the level of < +2% 
were present (see Sect. l3.1l for details). 



3. 1 . Luminosity-weighted single stellar population 

We compared our measured Lick indices Hj3, Mg b, Fe5270, 
Fe5335, Fe5105, Fe5406, and Fe5709 t o the val ues predicted by 
the single stellar population models of lThomas "et al. (2011). We 
interpolated the original grid of models so that we have steps 
of ~ 0.4 Gyr in age and ~ 0.025 dex in metallicity ([Z/H]) and 
in or-enhancement ([a/Fe]). The best stellar population param- 
eters are obtained by a minimization using the observed in- 
dices. Errors on the best-fit parameters are computed by means 
of Monte Carlo simulations. 

To quantify the errors on age, [Z/H] and [a/Fe] caused by 
residual sky contamination in the spectra (see Sect. 12.11 ) we 
added to the scientific spectra an extra contribution ±2% of the 
subtracted sky and measured the single stellar population pa- 
rameters again. The best-fit luminosity weighted values of age, 
[Z/H] and [a/Fe] and the range allowed by potential contamina- 
tion from sky emission lines are shown in Figure [3] and listed in 
Tabled 

There are no data in the literature that extend from the cen- 
ter of NGC 3311 toward the halo and overlap w ith our mea- 
surements, but according to iLoubser et al.l (l2008h . the central 
region of NGC 3311 is younger (Ageo= 8.7 +1.8 Gyr) con- 
siderably more metal-rich ([Z/H]o= 0.12 + 0.07) than what 
we observe in the hot halo, and has similar abundance ratio 
([Q'/Fe]o= 0.40 + 0.03). 

3.2. Contamination from foreground stars 

At ~ 30" north of the center of bin E there is a bright star with 
B ~ 12.2 mag from the USNO-Bl catalog but unknown spectral 
classification. Although this star does not enter the slit directly 
(its distance to the long-slit is ~20 times the seeing FWHM), it 
is inside the F0RS2 field of view and its light can be scattered 
along the optical path of the spectrograph and hit the detector 
In Sect.|2]we quantified the amount of scattered light (~ 7 ± 3% 
of the sky) from the portions of the detectors not covered by 
the long slit, and subtracted it from the scientific spectra. This 



ensures that the mean contribution from any scattered source is 
removed. The error in the quantification of the scattered light is 
very similar to the error in the sky subtraction, which does not 
explain the observed difference in stellar population properties 
(Fig. [3]). Small residual variations of the scattered light intensity 
along the long-slit might still be present, but it is not possible to 
measure them. Nevertheless, because the light is scattered, we 
expect that it contaminates adjacent regions of the detector by a 
similar amount. On the other hand, we observe a large difference 
between the spectral properties of the adjacent bins D and E, 
suggesting that this difference reflects a variation in the stellar 
population properties in different slit sections. 

4. Discussion 

NGC 3311 is surrounded by a dynamically hot stellar halo, 
character ized by a rapid inc rease of t he velocity disper sion of 
the stars ('Hau_gt al. 2004; Loubser_et al. 2008; V entimiglia et al.l 
2010) and of the globular clusters system (Ric htler et al. 1 1201 II) . 
In this section we will use the information on the stellar popula- 
tion derived in Sect. l3.1l to constrain the origin of this dynami- 
cally hot stellar halO; 

In this region. lArnaboldi et al.l (1201 ih identified a photomet- 
ric sub-structure (excess of light), which extends toward the 
northern and eastern sides of NGC 3311. In Fig. [T] we show the 
approximate contours of the excess of light in the field of view. 
An accurate determination of the morphology, photometry and 
precise isophotes of the excess is complicated by the presence 
of foreground stars and uncertainties in the background subtrac- 
tion and is not the aim of the present paper For our purposes 
it is enough to notice that the FORS2 long-slit crosses the ex- 
cess of light from SE to NW, probing regions where its intensity 
is different. Although the isophotes in Fig. [T] are approximate, 
we expect that the contribution from the excess of light to the 
mean stellar population properties along the line of sight has a 
minimum in bin A and a maximum in bin E. From Fig. [3] we 
have evidence that the mean stellar populations in bins E and D 
have different metallicities from those observed in the other slit 
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sections. We interpret this difference as caused by the contami- 
nation from the excess of light, which is higher in bin E (where 
the difference in metalUcity is more evident) and lower in bin D 
(where it is less evident). 

We therefore split the discussion into two parts: Section |4TT] 
describes the stellar population properties inferred from the por- 
tion of the slit covered by bins A-C, while Section l4~2l desci'ibes 
those derived from bins D-E. 

4.1. The dynamically hot stellar halo around NGC 331 1 

The spectra extracted in regions A-C contain the contribution 
mostly from stars in the dynamically hot halo around NGC 3311. 
Indeed, the rapid increase over radius of the velocity dispersion 
(see Section |4| suggests that the contribution of stars bound to 
the main body of NGC 331 1 at these galactocentric distances is 
small. From the contour levels displayed in Fig. [1] we can con- 
clude that the contribution from the excess of light is also small. 
Therefore, we consider the age, [Z/H] and [a/Fe] measured in 
bins A-C as indicative of the mean luminosity-weighted stellar 
population properties in the dynamically hot halo around NGC 
331 1, dominated by intra-cluster stars. 

In this region, stars are old (age > 13 Gyr), metal-poor 
([Z/H]= -0.35 + 0.05) and formed over a rapid time scale (< 1 
Gyr), as inferred from the high abundance ratio of a elements 
([a/Fe] - 0.48 + 0.03). Similar values are measured in the outer 
regions of ot her BCGs, e.g. NGC 3557 (Brough et al. 2007), 
NGC 1400 (Spolaoretal. '2008), NGC 1407'(SpolaoLeLalJ 
120081; iForbes et al. 2011). But, contrarily to NGC 3311, the ve- 
locity dispersion profiles in the halos of these galaxies are either 
decreasing or nearly constant, indicating that the stars are still 
bound to the galaxy and that the contribution from intra-cluster 
stars to the kinematics is small. The stellar halo around NGC 
331 1, however, consists of mostly intracluster stars. The lack of 
a significant metallicity gradient over the spatial range probed by 
slit sections A-C indicates that the stars accumulated on the halo 
through several accretion episodes. What are the progenitors of 
these stars? 

1 . Stars from the central regions of massive early type galaxies 
are most likely ruled out. The distribution of their [Z/H] ([a/Fe]) 
has mean 0.35 (0.19) and dispersion 0.30 (0.07) in solar units. 
These average values are based on measurements in early-type 
galaxies in th e Coma cluster wit h velocity dispersion cr > 200 
km s \ from lTrager et alj ([2008). The stellar population prop- 
erties measured in the central regions of the early-type galaxies 
in Hydra I where measurements of the Lick indices are avail- 
able QC 2586, IC 2597, NGC 3308, ESO436-G04, and ESO501- 
G13j3 are consistent with the average values in Coma. 

2. Stars in early-type galaxies outside 1-2 effective radii are po- 
tential candidates^Indeed, s everal ear ly-type galaxies in the sam- 
ple of lSpolaor et a l. (2010), PuetaDfeolO) have at 1-2 effective 
radii [Z/H] and [a/Fe] similar to what we observe in the halo 
of NGC 3311. Unfortunately, no data outside of one effective 
radius are available for early-type galaxies in Hydra I, but only 
central values for a number of them (NGC 3305, NGC 3309, 
ES0436-G45 and ESO501-G13pl. They are characterized by al- 
pha enhanced stars ([a/Fe]> 0.4) at their centers, but they have 
[Z/H]> 0. Stars stripped from their halos might be consistent 
with the stellar population of the stars around NGC 331 1 under 
the (plausible) assumption that these galaxies are characterized 

" See lJorgensenlll997l : lTrager et al.lll998l : IOgando et alj|2008l for the 
line strength of the Lick indices. 



by strong and negative metallicity gradients, and flat [a/Fe] ra- 
dial profiles. 

3. Stars stripped from the disks or the bulges of late-type spirals 
are most likely ruled out. They are generally characterized by 
younger ages and by solar values of metallicity and a-elements 
abundance (Crowl & Kenney 2008; Morelli et al. 2008). The 
stellar population properties measured in the central regions of 
the late-type galaxies in Hydra I where measurements of the 
Lick indices are available (ESO437-G015, ES0437-G45, and 
ESO501-G03pl are consistent with these values. A small frac- 
tion in mass of young stars would decrease the luminosity- 
weighted age of the corresponding single stellar population 
(ISerra & Tragedl2007h . in contrast with our measurements. 

4. Stars from disrupted dwarf galaxies are potential candi- 
dates. They are characterized on average by subsolar metal- 
licity and nearly solar a-enhancement; e.g., in the sample of 
Michielsen et al. (2008), the distribution of [Z/H] ([a/Fe]) in 
dwarfs has mean value -0.4 (0.0) with a scatter of 0.3 (0.2) 
dex. Therefore, stars from disrupted dwarf galaxies can lower 
the luminosity weighted metallicity of the equivalent single stel- 
lar population. On the other hand, their contribution would also 
decrease the luminosity- weighted [a/Fe], meaning that they can- 
not represent the main bulk of stars in the NGC 3311 halo. 

In conclusion, the combination of kinematic and stellar pop- 
ulation data suggests that the hot stellar halo around NGC 331 1 
was accreted. The bulk of the stars most likely originated from 
the outer regions of early-type galaxies with metal poor and a- 
enhanced stellar populations. Stars from disrupted dwarf galax- 
ies are also potential contributors, because of their generally low 
[Z/H]. 

4.2. The origin of the excess of light 

As discussed above, we have evidence that the stellar population 
measured in the long-slit bins D-E is different from that mea- 
sured in bins A-C. Following the same dynamical and photo- 
metric arguments as in Sect. 14.11 we consider the contribution 
of stars bound to NGC 331 1 to be small, but not the contribu- 
tion of the excess of light. Therefore, there are potentially two 
main stellar populations in this region, and the approximation 
of single stellar population adopted in Sect. 13. li ma v not be the 
appropriate way to interpret the observations. 

Despite of the S /N > 25 (per pixel) reached and the good 
quality of the spectra, we have been unable to separate the two 
populations by simultaneously fitting their kine matics and stel- 
lar populations (e.g., using a direct spectral fit as in lCoccato et al.l 
2011). This is probably because of the high velocity dispersion 
that broadens and mixes the spectral features of the two compo- 
nents along the line of sight. 

A better way to interpret our observations is therefore to 
make some assumptions on the kinematics of the excess of light 
and on the nature of its progenitors on the basis of available ad- 
ditional information. Following these assumptions, we created 
a simulated spectrum S model from the combined contribution of 
the stellar halo and the excess of light and tested whether it is 
consistent with what is observed. 

lArnaboldi et al.l (1201 lb were able to estimate the kinematics 
of the excess of light from the cross correlation between a tem- 
plate star and a model spectrum for the excess. They found two 
positive peaks on the coiTelation function above ~ 2cr noise. The 
first peak was at 5032 km s"' , and it was interpreted as the veloc- 
ity of the excess of light along the line of sight. The second was 
at 3905 km s"', consistent with the systemic velocity of NGC 
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33 11 and it was interpreted as the residual contamination of the 
stars from the halo in the model spectrum of the excess of light. 
Therefore the relative velocity of the excess of li ght with respect 
to NGC 3311 is on the order of 1200 km s '. lArnaboldi et all 
(1201 1 ) also propose a scenario in which the excess of light orig- 
inates from disrupted dwarf galaxies and they estimate that the 
excess of light represents 10% -20% of the surface brightness of 
the halo of NGC 33 1 1 in the region where the excess is detected. 

Following their interpretation, we assume that we observe 
the superimposition of the halo stellar population and a dwarf- 
galaxy-like population in bins D-E. As template for the halo 
spectra, we used that extracted from bin A; as template for 
a dwarf galaxy spectrum we used that of the dwarf galaxy 
HCC 26, on which the FORS2 long-slit was centered (see 
the orange circle in Fig. [TJ. The stellar population properties 
of HCC26 are consiste nt with those of other dwarf galaxies 
jMichielsen et al.|[2008l) . As a first approximation, we normal- 
ized the spectra of bin A and HCC 26 to their continuum 
values at 5100 A. We indicate these normalized spectra with 
Shiio and 5 dwarf, respectively. We then convolved 5 dwarf with 
a Gaussian LOSVD to construct the contaminant component 
5 dwarf- The LOSVD we used for the convolution has mean ve- 
locity 1200 km s"' higher than the systemic velocity of NGC 
3311 and velocity dispersion of 450 km s"' . These assumptions 
are consistent with the estimate of the velocity of the excess of 
light with respect the halo around NGC 3311 (lArnaboldi et al.l 
l201lh. the typical veloc ity dispersion measured in these regions 
dVentimiglia et al.ll2010 ) and the standard deviation of the line- 
of-sight velocities of the seven dwarf galaxies loc ated on the top 
of the excess of Ught (0-7 dwarfs ~ 470 km s~'. iMisgeld et al.l 
l2008h . which are potential progenitors We then constructed 
many combined spectra 5 model - (I - F) ■ 5 halo + F ■ 5 dwarf 
for diff'erent values of the "contaminating" fraction F, measured 
their equivalent single stellar population values of age, [Z/H] and 
[a/Fe], and compared them with the observed stellar population 
properties in bins D and E. We prefer this strategy over a di- 
rect comparison of the observed and combined spectra 5 model 
because it is less sensitive to the uncertainties of the overall 
shape of the stellar continuum. Lick indices are less affected by 
these uncertainties, because the pseudo-continua and the pseudo- 
bands where each index is defined are very close in wavelength. 

Figure |4] shows the results as function of F. The mean val- 
ues of the single stellar population parameters measured in bins 
D-E are consistent with those of a composite population made 
of stars from the stellar halo around NGC 3311 and stars from 
disrupted dwarf galaxies. In bin E the contamination is stronger 
than in bin D. In detail, the range of contamination in the surface 
brightness from 17% to 28% reproduces within uncertainties the 
mean values of age, [Z/H] and [c/Fe] measured in bin E (be- 
tween horizontal green long-dashed lines). A smaller fraction (< 
13%) is needed to explain the mean values observed in bin D 
(between horizontal magenta short-dashed lines). The ranges of 
contamination in surface brightness we found are on the sa me 
order of magnitude of the estimate in lArnaboldi et al.l (1201 ll) . 

The presence of a contaminant component with lower ve- 
locity dispersion (cr < 60 km s ') was also tested, and we found 
that it can contribute less than F <5% to the local surface bright- 
ness in all bins. Higher fractions would produce absorption lines 
features clearly visible on the combined spectrum (contrarily to 



' The results do not vary too much over the velocity dispersion range 
200 km s"'< cr < 600 km s"', meaning that we cannot infer precise 
constraints on the kinematics using the stellar population properties. 
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Fig. 4. Equivalent single stellar population age (upper panel), 
metallicity (central panel) and a-enhancement (lower panel) 
measured in 5 model- model represent a composite stellar popu- 
lation formed by stars from the halo around NGC 331 1 contam- 
inated by a fraction F of stars from dwarf galaxies (see Sect. 14.21 
for details). Long-dashed green and short-dashed magenta lines 
represent the mean values of age, [Z/H] and [o'/Fe] measured 
in the long-slit bins E and D, respectively. The range of F val- 
ues consistent with measurements in bins E and D are shaded in 
green and magenta, respectively. 



what we observe in bins D-E), which would be recovered by a 
direct spectral fit dCoccato et alj|20l"l1) . 

Despite of the uncertainties on the kinematics that we as- 
sumed for the contaminant stellar component, this simple ap- 
proach supports the idea that the excess of light in the Hydra I 
cluster core consists of stars stripped from dwarf galaxies, which 
are thus being added to the dynamically hot stellar halo around 
NGC 3311. 



5. Summary 

We used deep long-slit spectra obtained with FORS2 at the VLT 
on the northeast side of the brightest cluster galaxy NGC 33 1 1 to 
study the stellar population content in five regions from ~10 to 
~30 kpc from its center (slit sections A-E), where the co ntribu- 
tion from intra-cluster stars dominates dVentimiglia et al.ll2010l) . 

We found that the stars of the dynamically hot halo around 
NGC 3311 are on average characterized by an old, metal-poor 
and a-enhanced stellar population. The lack of a stellar popu- 
lation gradient in slit sections A-C together with the kinematics 
suggest that the intracluster halo formed from accretion of mul- 
tiple components onto the cluster center The stellar population 
parameters suggest that the halo stars were mainly stars stripped 
from the outer regions of early-type galaxies, with a possible 
contribution from dwarf galaxies. 

In the long-slit sections D-E the metallicity is significantly 
lower than that measured in sections A-C. We associate this 
spectroscopic feature with the presence of a photometric sub- 
structure (excess of light), which wa s detected in the sam e 
area of the halo around NGC 3311 dArnaboldi et all 1201 ih . 
Simulations of the measurements over this region are consis- 
tent with the presence of a composite population made of stars 
from the stellar halo around NGC 331 1 and stars from disrupted 
dwarf galaxies, which form part of the excess of light. Our find- 
ings therefore support the idea that the build-up of ICL in the 
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Hydra I cluster core and the stellar halo around NGC 3311 are 
still ongoing. 
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